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On-line  measurements  of  temperature  and  concentration 
and  the  global  rate  at  any  reference  state  of  catalyst  can  be 
used  to  estimate  a  measure  of  catalyst  deactivation.   A  feed- 
back control  policy  is  obtained  from  this  measure  of  cata- 
lyst deactivation  for  manipulating  the  reactor  outlet  con- 
version or  yield  in  any  desired  manner.   This  feedback  con- 
trol enables  one  to  obtain  the  desired  conversion  or  yield 
without  any  knowledge  of  deactivation  kinetics,  and  the  re- 
sults are  applicable  to  any  type  of  deactivation.   The  feed- 
back control  policy  has  been  verified  by  experiment  based 
on  ethylene  oxidation  reactions  over  a  silver  catalyst  using 
a  laboratory  shell-and-tube  reactor.   Excellent  agreements 
have  been  reached  between  the  theory  and  the  experiment. 


The  problem  of  finding  the  optimal  piecewise  control 
policy  has  been  formulated  as  a  multivariable  optimization 
problem.   The  optimal  policy  is  one  of  increasing  conver- 
sion.  The  constant  conversion  policy  corresponds  to  the 
limiting  case  where  the  number  of  piecewise  steps  approaches 
infinity.   This  optimal  piecewise  control  is  suitable  for 
on-line  implementation  in  conjunction  with  the  on-line  es- 
timation of  the  extent  of  catalyst  deactivation.   An  optimal 
continuous  control  policy  for  parallel  and  reversible  reac- 
tions has  also  been  solved  using  the  technique  of  calculus 
of  variations.  It  is  shown  that  a  direct  substitution  method 
is  more  efficient  than  the  commonly  used  method  of  applying  the 
Pontryagin  maximum  principle. 

The  optimal  control  policies  can  easily  be  taken  into 
consideration  for  reactor  design.   This  combination  of  pro- 
cess design  and  control  with  due  consideration  of  catalyst 
deactivation  for  both  is  shown  to  result  in  a  substantial  im- 
provement of  the  reactor  performance. 
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CHAPTER  1 
INTRODUCTION 

The  usual  practice  of  operating  a  fixed-bed  with  slow 
catalyst  deactivation  is  to  raise  the  reactor  inlet  tem- 
perature to  compensate  for  the  declining  catalytic  activi- 
ty so  as  to  maintain  a  desired  conversion  (Kovarik  and  Butt, 
1982)  .   The  catalyst  is  regenerated  when  the  required  tem- 
perature becomes  too  high.   The  reactor  is  usually  designed 
to  give   the  desired  conversion  without  any  consideration 
for  the  catalyst  deactivation.   A  summary  of  operational 
strategies  for  a  reactor  subject  to  catalyst  deactivation 
is  listed  in  Table  1-1. 

A  number  of  questions  can  be  raised  regarding  these 
usual  practices:   The  first  question  is  how  we  should  con- 
trol the  operating  conditions  of  a  reactor  designed  in 
the  usual  way.   This  question  has  been  the  subject  of  many 
studies  (Chou  et  al.  1967;  Szepe  and  Levenspiel,  1968;  Ogunye 
and  Ray,  1971;  Haas  et  al.,  1974;  Levenspiel  and  Sadana, 
1978) .   The  control  policies  resulting  from  these  studies 
either  lead  to  an  open-loop  control  or  require  detailed 
knowledge  of  catalyst  deactivation.   Although  considerable 
progress  has  been  made  in  our  understanding  of  catalyst 
deactivation,  the  deactivation  is  still  the  least  understood 
of  all  facets  involved  in  the  quantitative  description  of  a 


Table  1-1   Summary  of  Operational  Strategies  for  a  Reactor 
Subject  to  Catalyst  Deactivation 


1.  Vary  reactor  temperature  with  time  to  maintain  a  constant 
conversion  with  a  constant  reactor  feed  flow  rate.   A  typi- 
cal policy  for  large  throughput  (400  MM  pounds  per  year)  and 
slow  deactivation  rates  (months  to  years  of  catalyst  life) . 

2.  Vary  throughput  of  the  reactor  feed  while  holding  the  reactor 
temperature  and  conversion  constant.   A  possible  policy  for 
medium  deactivation  rates  (weeks  to  months  catalyst  life) 

and  small  to  medium  throughput  (about  25  MM  pounds  per  year) 
systems. 

3.  Allow  the  conversion  to  fall  while  holding  the  reactor  feed 
flow  rate  and  reactor  temperature  constant.  Similar  appli- 
cations as  in  Item  2 . 

4.  Maintain  the  fresh  feed  rate  and  reactor  temperature  constant 
and  let  the  recycle  flow  increase.   Similar  application  as 

in  Item  2 . 

5.  Use  a  combination  of  reactors  in  parallel  and  the  policies  of 
Items  1  or  3.   Usually,  with  two  reactors  in  parallel,  one 
will  be  off-line  for  catalyst  regeneration  while  the  other  is 
operating.   A  typical  policy  for  large  throughput  and  medium 
to  fast  deactivation  rates  (days  to  months  of  catalyst  life) . 

6.  Continuous  catalyst  regeneration  while  maintaining  constant 
conversion,  throughput,  and  reactor  temperature.  A  typical 
policy  for  large  throughput,  rapid  deactivation  systems 

(hours  to  days  of  catalyst  life)  . 

Reprinted  from  Cat.  Rev.  .24,  499  (1982),  by  courtesy  of  Marcel 
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fixed-bed.   This  uncertainty  regarding  the  deactivation  and 
the  compounded  effects  of  deactivation  and  diffusion  in  a 
fixed-bed  make  it  quite  unattractive  to  implement  an  open- 
loop  control  policy.   The  second  question,  therefore,  is 
whether  a  measure  of  catalyst  deactivation  can  be  estimated 
from  process  measurements  and,  if  so,  how  a  feedback  control 
can  be  realized  with  the  measurements  made.   In  Chapter  2 
a  method  of  determining  a  measure  of  catalyst  deactivation 
is  presented  and  the  result  then  utilized  for  a  feedback 
control  policy.   This  feedback  control  policy  is  demonstrated 
in  Chapter  3  with  a  simulation  of  a  model  reaction  system, 
and  is  put  to  the  test  with  a  laboratory  fixed-bed  where 
ethylene  oxidation  reactions  take  place  over  a  silver  cata- 
lyst. 

It  has  been  shown  that  constant  conversion  policy  is 
optimal  only  for  a  single  irreversible  reaction  affected     , 
by  concentration-independent  deactivation,  when  the  reactor 
is  controlled  continuously  (Chou  et  al.  1967;  Lee  and  Crowe, 
1970;  Crowe,  1970).   Thus,  the  third  question  regarding  the 
usual  practice  of  maintaining  constant  conversion  is  what 
the  optimal  policy  will  be  if  the  reactor  is  piecewise  con- 
trolled (for  example,  the  temperature  is  adjusted  only  once 
a  day,  rather  than  continuously),  or  if  the  reaction  is 
complex.   The  solutions  to  this  question  are  presented  in 
Chapters  4  and  5. 


A  much  more  important  question  than  the  three  discussed 
above  has  to  do  with  the  inherent  interrelationship  be- 
tween process  controllability  and  process  design.   It  is 
intuitively  clear  that  the  process  design  dictates  the 
controllability,  for  the  parameters  involved  in  describing 
a  process  contain  design  parameters  which  in  turn  dictate 
the  way  the  process  can  be  controlled.   It  is  clear 
therefore,  that  the  best  possible  performance  of  the  process 
can  be  attained  when  the  process  design  and  control  are 
combined.   For  the  reactor  under  consideration,  this  means 
that  the  best  possible  performance  can  be  attained  when  the 
catalyst  deactivation  is  taken  into  consideration  not  only 
for  the  reactor  control  but  also  for  the  reactor  design. 
The  question,  therefore,  is  what  the  reactor  size  and  the 
feedback  control  should  be  for  the  best  performance  of  the 
reactor.   In  Chapter  6,  the  answer  to  this  question  is  pre- 
sented. 


CHAPTER  2 
THEORY  OF  ON-LINE  ESTIMATION  OF  DEACTIVATION  AND  CONTROL 

An  important  problem  in  operating  a  fixed-bed  with 
catalyst  deactivation  is  that  of  finding  a  means  of  esti- 
mating the  extent  of  catalyst  deactivation  from  on-line 
measurements  because  of  the  uncertainty  regarding  the  deac- 
tivation.  The  problem  of  controlling  the  conversion  in  any 
desired  manner  using  a  feedback  scheme  is  another.   These 
two  problems  are  solved  theoretically  in  this  chapter, 
starting  with  a  single  reaction  and  then  proceeding  to 
multiple  reactions. 
§2-1   A  Simple  Reaction 

Consider  a  shell-and-tube  reactor  in  which  a  single 
reaction  takes  place.   Assuming  negligible  radial  gradients 
and  pseudo  steady-state,  one-dimensional  balance  equations 
can  be  written  as 


f  =  -TV   T  -  ZU-eB)/v 
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where  C  and  T  are  the  concentration  and  temperature  of  the 
reactant,  and  T   is  the  countercurrent  coolant  temperature. 
Here  m  and  m  are  the  reaction  fluid  and  coolant  mass  rates, 
C  and  (C  )   are  the  specific  heat  capacities  for  the  re- 
action  fluid  and  the  coolant,  Z  is  the  reactor  length,  z 
is  the  axial  reactor  coordinate  normalized  with  respect  to 
Z,  v  is  the  superficial  velocity  assumed  constant  in  the 
radial  direction,  eB  is  the  bed  porosity,  U  is  the  overall 
heat  transfer  coefficient,  and  R  is  the  tube  radius. 
The  global  rate  based  on  pellet  volume,  R_ ,  can  be  expressed 
as 

Rr  =  h(z;t)k  f (C,K.)  (2-4) 

*J  a. 

where  k   is  the  apparent  rate  constant  in  the  form  of  the 

d 

Arrhenius  relationship  and  f  represents  the  apparent  depen- 
dence of  the  rate  on  concentration  and  equilibrium  constants 
K  .  .   In  the  absence  of  diffusion  and  deactivation  effects, 
the  apparent  rate  reduces  to  intrinsic  rate.   The  activity 
factor  h  (Wheeler,  1955)  is  defined  by 

h(z;t)  E  — GD  (2_5) 

1  G'r 

where  the  subscripts  r  and  D  denote  reference-state  and  de- 
activated catalyst,  respectively.  According  to  the  defini- 
tion, h  is  unity  for  catalyst  at  the  reference  state.  This 
activity  factor  decreases  with  time  due  to  catalyst  deacti- 
vation and  thus  depends  on  time.  It  should  be  recognized 
in  Eq.  (2-5)  that  the  definition  of  h  is  general  in  that  h  can 


depend  on  the  concentration  C  and  thus  is  applicable  to  con- 
centration-dependent deactivation.   Since  C  is  a  function  of 
z  at  a  given  time,  h(C;t)  is  a  subset  of  h(z;t).   The 
pseudo  steady-state  assumption  is  based  on  the  fact  that 
the  time  scale  of  catalyst  deactivation  is  many  orders  of 
magnitude  larger  than  that  for  the  reactor  to  reach  a  steady 
state. 

For  the  purpose  of  deriving  an  expression  for  a  measure 
of  the  extent  of  deactivation,  we  rewrite  Eq.  (2-1)  with 
the  aid  of  Eq.  (2-4)  as 

||  =  -Th(z;t)kaf  (C,K.)  (2-6) 

Here  again,  one  can  write  h  as  h(C;t)  in  place  of  h(z;t) 

for  dependent  deactivation.   This  equation  can  be  integrated 

from  the  reactor  inlet  to  the  outlet  to  give 
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h(z;t)dz  =  H(t) 
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(2-7) 


where  Ce  and  Cin  are  the  concentrations  at  the  outlet  and 
the  inlet,  respectively,  at  any  given  time.   By  definition, 
H  is  unity  at  reference  state  since  the  activity  factor 
h  is  uniform  at  unity  throughout  the  reactor.   The  quantity 
H  is  a  measure  of  the  extent  of  catalyst  deactivation  which 
represents  the  activity  for  the  reactor  and  thus  may  be  termed 
"reactor  activity  factor."   The  value  of  H  decreases  as  the 
outlet  conversion  decreases  due  to  the  catalyst  deactivation 
for  a  given  inlet  temperature. 
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The  temperature  appearing  in  Eq.  (2-7)  needs  to  be  re- 
lated to  concentration  for  the  integration.   Integrating 
Eq.  (2-2)  from  the  reactor  inlet  to  z,  we  obtain 
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in 
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where  T is  the  coolant  temperature  at  z=0  and  T.   is  the 

ce  c  m 

reactor  inlet  temperature.   Consider  first  an  adiabatic  re- 
actor for  which  Eq.  (2-8)  reduces  to 
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in 
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Now  that  T  is  expressed  in  terms  of  C  and  readily  measur- 
able inlet  conditions,  the  value  of  H  can  be  calculated  by 
simply  carrying  out  the  integration  numerically  with  respect 
to  C  in  Eq.  (2-7)  with  the  measured  outlet  concentration. 
For  instance,  the  apparent  rate  constant  k  (T)  for  an  adia- 
batic  reactor  can  be  expressed  as 


k  =  k   exp 
a    a    * 
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where  k    is  the  apparent  preexponential  factor  and  E   is  the 

o.  a 

apparent  activation  energy.   The  equilibrium  constants  K. 

can  be  expressed  in  a  similar  manner.   For  a  shell-and-tube 
reactor  to  which  Eq.  (2-8)  applies,  the  calculation  of  H  re- 
quires measurement  of  coolant  temperature  along  the  reactor 
length.   This  measurement  can  be  used  in  Eq.  (2-3)  for  the 


calculation  of  T  as  a  function  of  z.   With  these  T  and  T  , 
the  corresponding  change  in  C  for  a  selected  interval  of  Az 
can  be  obtained  from  Eq.  (2-8)  which  in  turn  can  be  used  in 
Eq.  (2-7)  for  the  calculation  of  H. 

An  immediate  use  of  the  reactor  activity  factor  calcul- 
able from  the  temperature  and  concentration  measurements  is 
in  the  manipulation  of  the  reactor  inlet  temperature  for  the 
purpose  of  maintaining  the  desired  conversion.   Suppose  that 
the  inlet  temperature  is  adjusted  intermittently,  as  in  the 
usual  operation  of  a  reactor,  to  compensate  for  the  declining 
activity.   Consider  an  adiabatic  reactor  for  clarity.   If  we 
let  the  subscript  c  denote  the  current  quantities  and  n  the 
new  quantities  resulting  from  a  change  in  the  reactor   inlet 
temperature,  Eq.  (2-7)  can  be  written  twice  to  give 


H   = 
c 


H   = 
n 


dC 


h  (z;t)dz  -  -  .         [T;(T.  )  ]f(C,K.) 
0  Cin  in  c       3 


1 

h  (z;t)dz  =  - 

0  n 


Cd  dC 


C   Tka[T?(Tin)nlf(C'Kj) 

in 


(2-11) 


(2-12) 


where  C,  is  the  desired  concentration  which  will  be  attained 

by  changing  the  current  inlet  temperature,  (T^n) c,  to  the 

new  inlet  temperature,  (T.  )  .   As  we  shall  soon  see   in 

the  next  section,  the  activity  factors  h   and  h  depend  main- 

c      n 

ly  on  the  fraction  of  catalyst  deactivated  and  slightly  on 
the  temperature  difference  between  bulk  fluid  and  pellet 
surface.   If  we  neglect  this  change  in  the  temperature 
difference  for  the  time  being,  we  can  set 

hc(z;t)  -  hn(z;t)  (2-13) 
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for  the  purpose  of  calculating  the  new  inlet  temperature 
that  yields  the  desired  outlet  concentration.   The  justifi- 
cation here  is  that  the  fraction  of  catalyst  deactivated 
changes  negligibly  for  all  points  along  the  reactor  while 
the  inlet  temperature  is  adjusted  from  the  current  value  to 
a  new  value.   Note  that  the  justification  is  still  valid  for 
concentration  dependent  deactivation,  i.e.   h  (C(z);t) 
=  hn(C(z);t).   Equations  (2-11)  through  (2-13)  together  with 
Eq.  (2-10)  can  be  used  to  obtain 


fC, 


dC 


Cin   Tka  ex? 
o 


-E  /R 
§. 2. 


(T.  )  + 
in  n 


-AH 


PC 

1      PJ 


(c,  -  C) 


in 


=  H 


f(C,K.) 


(2-14) 
where  Hc  is  the  known,  current  reactor  activity  factor  cal- 
culated from  Eq.  (2-12)  just  prior  to  the  change  in  the  in- 
let temperature  based  on  the  current  measured  outlet  con- 
centration Ce  and  inlet  temperature  (T.  )  .   Equation  (2-14) 
can  now  be  solved  for  (T.  )   by  a  numerical  technique  such 
as  Newton's  method  or  bisection  method  since  the  left  hand 
side  of  the  equation  can  be  evaluated,  given  a  value  of 
^inV   The  value  of  <Tin>n  that  satisfies  Eq.  (2-14)  will 
yield  the  desired  conversion  and,  therefore,  the  inlet  tem- 
perature can  be  changed  to  (Tin)   to  maintain  the  conversion 
at  the  desired  level. 
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The  control  policy  of  Eq.  (2-14)  is  a  piecewise  algo- 
rithm for  the  manipulation  of  the  inlet  temperature.   There- 
fore, it  will  maintain  the  conversion  at  the  desired  level 
only  for  a  short  period  of  time  after  a  change  in  the  inlet 
temperature  and  the  conversion  thereafter  will  gradually  de- 
crease with  time  until  the  inlet  temperature  is  raised  again. 
A  bandwidth  for  the  allowed  decrease  in  conversion  can  be  used 
to  trigger  the  adjustment  of  the  reactor  inlet  temperature. 
The  control  algorithm  for  a  nonadiabatic  reactor  can  be  ob- 
tained in  a  similar  manner.   The  equation  corresponding  to 
Eq.  (2-14)  is 

^d  dC 

"       Tk  [T;(T.  )  ]f(C,K.)  "  Hc  (2_15) 

cin   a     in  n      3 

This  equation,  however,  needs  to  be  solved  in  conjunction 
with  Eqs.  (2-3)  and  (2-10)  for  (T^  )   along  with  the  measured 
Tc  profile. 

As  indicated  earlier,  we  assumed  that  the  heat  transfer 
resistance  across  the  interface  between  bulk  fluid  and  pellet 
surface  is  negligible  in  arriving  at  the  control  policies  of 
Eqs.  (2-14)  and  (2-15).   If  there  exists  a  significant  tem- 
perature difference  between  bulk  fluid  and  pellet  surface, 
Eqs.  (2-13)  through  (2-15)  no  longer  hold.   The  discrepancy 
between  hc  and  hn  caused  by  the  temperature  difference  can 

be  compensated  for  by  adjusting  C,  as  follows: 

d 
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Cd.  =  Cd  +  a(Cd    ~  Ce    >      i-1,2,...  (2-16) 

i  l-l    i-1 


do    d 


where  Cd   is  the  adjusted  C  when  a  step  change  in  the  in- 

i 
let  temperature,  which  is  triggered  by  the  bandwidth  con- 
straint, is  made  at  the  ith  step  and  a  is  a  proportionality 
constant  with  a  value  between  zero  and  unity.   Thus,  for 
a  given  decrease  in  the  outlet  concentration  from  the  de- 
sired value,  the  piecewise  feedback  control  algorithm  can 
be  stated  as  follows: 

£d. 

1  dC 

C    Tk  [T;(T   )  ]f  (C,K.)   =  IV  i-1 
in    a      in  1       3 


i-1, 2, . 


(2-17) 


where  (H<2)i_1  is  calculated  from 


(H  )  .  , 

c  1-1 


c  o 


C 

r  e.  , 

i-l  dC 


C. 
in 


k  [T; (T.  ) .   if (C,K.)  i-1/2/-- •  (2-18) 
a      in  i-l       -i 


with  the  measured  outlet  concentration  and  inlet  tempera- 
ture.  Here,  Cd   is  given  by  Eq.  (2-16).   Each  time  the  outlet 

i 
concentration  reaches  the  allowed  bandwidth,  the  new  inlet 
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temperature   (T.  ) .  is  calculated  for  the  manipulation  of 
the  inlet  temperature.   All  that  is  required  for  the  feed- 
back control  is  the  global  rate  for  catalyst  at  a  reference 
state. 

§2-2   Nature   of  Reactor  Activity  Factor 

Before  we  proceed  to  multiple  reactions,  let  us  examine 
the  nature  of  the  reactor  activity  factor.   By  definition, 
it  is  an  integrated  value  of  local  activity  factor  (Eq.  2-7), 
The  local  activity  factor,  which  is  the  activity  factor  at 
a  point  in  the  reactor,  can  be  obtained  by  solving  pellet 
conservation  equations.   Instead  of  solving  the  conservation 
equations,  we  utilize  the  global  rate  obtained  by  Lee  and 
Butt  (1982)  for  a  reaction  affected  by  uniform  deactivation 
and  diffusion  to  get  a  clear  picture  of  the  local  activity 
factor.   The  global  rate  for  this  case  is 


ou„  -    — <^ 


G  D    L       _  _  (2-19) 

[2D  kl(C)  J3*    <* 


where 


I-  = 


rC 

g(a)da  (2-20) 

0 


1.2E  (-AH) [2D  kl]3* 

q  =  5 — ~ <2-21> 

2h  R  T 
m  g 
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Here,  y   is  the  fraction  of  catalyst  deactivated,  g(C)  the 
concentration  dependence  of  the  intrinsic  rate  of  reaction, 
k  the  intrinsic  rate  constant,  De  the  effective  diffusivity, 
L  the  characteristic  pellet  length,  and  hm  the  film  heat 
transfer  coefficient.   The  local  activity  factor  can  be  ob- 
tained from  the  definition  of  Eq.  (2-5) : 


h  = 


(Vr 


<Vd 


<Vd 


Y=l 


=  (1-Y) 


Iza. 


_1-(1-Y)  2q 


(2-22) 


The  same  results  hold  (Lee  and  Ruckenstein,  1983)  for  non- 
uniform deactivation  if  Y  is  replaced  with  y    ,  which  is  y 
at  the  pellet  surface,  as  long  as  y   is  less  than,  say,  0.5. 
It  is  clear  from  Eq.  (2-22)  that  the  local  activity  factor 
depends  on  Y  and  the  quantity  q  which  represents  the  heat 
transfer  resistance  across  the  bulk  fluid-pellet  surface  in- 
terface.  Since  y  changes  very  little  while  the  inlet  tem- 
perature is  raised  to  a  new  value,  the  magnitude  of  q  de- 
termines whether  hc  can  be  set  to  hR.   Therefore,  the  pro- 
portionality constant  a  in  Eq.  (2-16)  can  be  chosen  in 
accordance  with  the  magnitude  of  q:   for  q  much  smaller 
than  unity,  the  value  of  a  can  be  set  to  zero,  its  value  in- 
creasing with  increasing  q.   Although  the  conclusion  is  made 
from  the  result  for  independent  deactivation,  the  same  should 
be  valid  for  dependent  deactivation  since  the  factor  causing 
the  discrepancy  between  hQ  and  hR  is  still  the  same,  i.e. 
the  heat  transfer  resistance  across  the  interface,  which  q 
represents. 
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§2-3  Multiple  Reactions 

For  simple  parallel  reactions,  we  can  write 


b  B 


A  -  b2C 

where  b1  and  b2  are  the  ratios  of  stoichiometric  coeffi- 
cients.  We  restrict  our  attention  here  to  an  adiabatic 
reactor  for  brevity  since  a  nonadiabatic  reactor  can  be 
treated  in  the  same  manner  as  shown  for  a  single  reaction, 
Let  6  denote  the  ratio  of  local  activity  factors,  i.e., 
8=n2 (z; t) /h^ (z;t) ,  the  reactor  conservation  equations 
can  be  written  as: 

dCA 

dz"  =  "  TV«'t,{ka   fl<CA'  CB>  +  *  ka_  f2(CA'CB)}; 


=  C, 


z=0 


in 


(2-23) 


dC 


B 


dz 


=  b1xh1  (z;t)k    f.  (C_,  CD)   ;   C 
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ax   lv  A'   BJ 


B 


=  C        (2-24) 

z=0    ain 


T  =  Tin  +  al<CA,.  "  V  +a2(CB  "  CB.  > 


in 


in 


(2-25) 


where  a1  =  -  AH2/pCp  and  a2  =  -  AH^  (l-A^/AR^)  pC  b..  Equa- 
tion (2-25)  results  when  the  heat  balance  equations  are  combined 
with  the  mass  balance  equations.   Here  the  global  rates  for 
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the  first  and  second  reactions  are  expressed  as 

hka  fl(CA'CB)  and  3hka  f2 ^Ch'CB^  '    resPectivelY-   If  the  same 
active  sites  are  responsible  for  both  reactions,  the  local 
activities  are  the  same,  and  8  is  unity.   Otherwise,  we 
assume  that  8  is  a  function  of  time  but  constant  along  the 
reactor,  i.e.,  3=3 (t). 

The  reactor  activity  factor  can  be  obtained  in  the  same 
manner  as  for  a  single  reaction: 
■1 


Hl  = 


rCB, 


h.dz  = 


dC 


B 


V^WV 


(2-26) 


'B 


in 


where  Cg   is  the  outlet  concentration  of  species  B.   As  Eq. 

e 
(2-25)  indicates,  the  temperature  can  be  expressed  in 

terms  of  CA  and  CB  for  the  calculation  of  H  but  C  has  to  be 

related  to  CB  to  carry  out  the  integration.   Combining  Eqs . 

(2-23)  and  (2-24),  we  have 


dC, 


dC 


B 


1_ 

b. 


1+8 


ka2f2(CA'CB> 
Vl^A'V 


;  c, 


=  c. 


c  =  c 

B    S, 


in 


in 


(2-27) 
which  can  be  solved  to  obtain  C   as  a  function  of  C  , 

A  B 

CA=  f(cB)'  bY  Fourth  Order  Runge  Kutta  Method  or  Predictor- 
Corrector  Method  (Hornbeck,  1975).   If  8  is  unity,  the  so- 
lution is  straightforward.   Otherwise,  the  concentrations 
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measured  at  the  outlet,  CA  and  C   ,  can  be  used  as  an 

e       e 


C   =  C   =  Ca  '  for 
B    UB     e 
e 


additional  boundary  condition,  i.e.,  C 

Eq.  (2-27)  to  obtain  the  value  of  6. 

If  species  B  is  the  desired  species  such  that  (C  )  n 

B  d 

is  the  desired  outlet  concentration,  the  control  policy 
similar  to  Eq.  (2-14)  can  be  written  as 


r(CB>d 


Hic= 


dC 


B 


B.     b, xk    exp 
in     la,     ^ 
lo 


-E  /R 
a  _g_ 


(Tin)n+  ai(CA.  "  CA)+a2(CB-CB. 
in  in 


fl(CA'CB> 


(2-28) 


which  follows  from  Eq.  (2-24).   Here  again,  H   is  the  re- 
actor activity  factor  for  the  reaction  path  1  calculated  just 
prior  to  a  change  in  the  inlet  temperature.   Equation  (2-28) 
can  be  solved  for  (Tin>n  with  the  aid  of  Eq.  (2-27).   It  is 
noted  that  (T±n)n  should  also  be  used  in  Eq.  (2-27)  in  the 
calculation. 

Consider  now  consecutive  reactions,  for  which  we  have 

1  2 


A  ►  b1   B   »■  b  C 

With  the  same  assumptions  made  for  the  parallel  reactions, 
the  conservation  equations  for  an  adiabatic  reactor  are 


5 

dz 


■Th(Z;t)ka1fl(CA'CB,;  CA 


=  C, 


z=o 


in 


(2-29) 
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dC 


B 


dz 


•Th(Z;t)  {eka2f2(cA,cB)  -v    f^C^)} 


'B 


=  C 


z=0 


B. 


in 


(2-30) 


T=Tin  +a3(CA.  "  CA>"  VCB  "  CB .  > 

in  m 


(2-31) 


where  a3  =  (-AH^-bAH^)  /pC   and  a4=-AH2/pC  .   If  we  follow 
the  same  procedures  as  for  the  parallel  reactions,  we  have 
for  the  reactor  activity  factor: 


Hl  = 


fC, 


h..dz  =  - 


dC, 


Tkai(T)fl(CA,CB) 


(2-32) 


in 


The  relationship  between  CA  and  CB  resulting  from  Eqs 
<2-29)  and  (2-30)  is 


dC 


dC 


-  -  -b,+  B 


ka1fl(CA'CB)  ' 


'B 


=  CB      (2-33) 

CA  =  CA.      in 
m 


The  same  procedures  as  for  the  parallel  reactions  can  be 
used  to  calculate  H  from  Eq,  (2-32)  with  the  aid  of  Eqs.  (2-31) 
and  (2-33). 

The  control  policy  for  maintaining  the  desired  out- 
let concentration  of  species  B  follows  directly  from  Eq. 
(2-30)  and  the  definition  of  H: 


Hic=- 
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"Vd  dcB 


c     n»ka21|Ii.,»lf2ICA'c.|-V1/"i,V¥C*C,') 

B  ■ 
in 

(2-34) 
This  can  be  solved  for  (T.  )   with  the  aid  of  Eqs .  (2-31) 
and  (2-33).   If,  however,  the  control  objective  is  to  have 
maximum  CB  at  the  outlet,  the  choice  of  (C  ),  in  Eq.  (2-34) 
must  satisfy  the  necessary  condition: 


dCB 


dz 


=  0  (2-35) 

z=l 


Substituting  Eq.  (2-35)  into  Eq.  (2-30)  gives 

ka2(Te)f2CCA'(CB>d]  =  ^Va,  (Te>  f  1  [CA  <  {CB^  <2"36> 

*         e  1         e    • 

which  can  be  used  to  obtain  (Cg),  with  the  aid  of  Eqs. 

(2-31)  and  (2-33).   The  new  inlet  temperature,  (T.  )    can 

in  n 

then  be  obtained  from  Eq.  (2-34). 

The  control  policies  of  Eqs.  (2-28)  and  (2-34)  can  be 
modified  using  Eq.  (2-16)  to  obtain  the  piecewise  algorithms 
similar  to  Eqs.  (2-17)  and  (2-18). 


CHAPTER  3 
SIMULATION  AND  EXPERIMENTAL  VERIFICATION 
The  theory  in  Chapter  2  is  demonstrated  here  with  a 
simulation  of  a  model  reaction  system.   The  theory  is  also 
put  to  the  test  with  a  laboratory  shell-and-tube  reactor  in 
which  ethylene  oxidation  reactions  take  place  over  a  silver 
catalyst.   The  catalyst  undergoes  sintering  under  the 
reaction  conditions. 
§3-1   Simulation  of  a  Model  Reaction  System 

For  an  illustration  of  the  calculation  of  H  and  its 
use  for  the  reactor  control,  we  consider  the  reaction  system 
studied  by  Lee  and  Butt  (1982a) ,  which  is  summarized 
in  Tables  3-1  and  3-2.   For  the  reaction  taking  place  in  an 
adiabatic  reactor,  which  is  affected  by  diffusion  and  uni- 
form, independent  poisoning,  the  intrinsic  rate  of  deacti- 
vation rp  is  given  (Lee  and  Butt,  1982b)  by 


rp  =  kp  (1_^N 


where  kp   is  the  rate  constant  evaluated  at  the  pellet  sur- 
face  temperature  and  N  is  the  concentration  of  poisoning 
species.   With  the  use  of  reactor  point  effectiveness  (Lee 
and  Butt,  1982a) ,  the  reactor  behavior  can  be  simulated  in 
a  straightforward  manner.   The  simulation  results  are  used 
here  as  the  process  output  for  the  purpose  of  illustrating 
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Table  3-1   Equations  for  the  Model  Reaction  System 


dC 
dz 


=  -   TR, 


a  =  C. 
z=0     in 


dN 
dz 


-   TkpaN(l  -  y)  ; 


N 


z=0  :  Nin 


dy_     1 

dt  "   Q  rp 


"  Q  kpaNU  "  Y);   Y 


t=0 


=  0 


where  R_  = 


1/L 


1.2  E  (-AH) 


2Dek(l-Y) 


fC 


^h 


g(c)dc 


2h  R  T' 
m  g 


g(c)  = 


1+CC 


T  =  T. 
in 


-AH 


IPC 


(C  -  C.  ) 

in 


k  =  k  exp(-E  /R  T) 
o    c        a   g 


T   =  T  +  (-AH)  RrL/hm 
s  G   m 


k    =  k exp  (-E  /R  T  ) 

pa    po        p  g  s 
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Table  3-2   Parameters  For  The  Model  Reaction  System 


k  =  exp 


12,000 


+  14.6 


(1/s);  E^  =  23.76  kcal/mol 

3. 


K  =   exp 


3,600 


+  3.86 


(cm  /mol) 


Cin  =  1.81  x  10  5  mol/cm3 


f-AH 


PC 


7      3 
=  4  x  10   K  cm  /mol;   -AH  =5.04  kcal/mol 


hm  =  3.245  x  10~4  cal/s.cm2.K 


N±n  ■  5  x" 10~8  mol/cm3 

-4        3 
Q  =  10   mol/cm  cat.  pellet' 


kp  =  exp 


7,000 


+  3.27 


;  E  =  13.86  kcal/mol 

IT 


D   =  10~3  cm2/s 

ir 


Tmav   =  1280K 
max 


X  =  70%     AX  =  10! 


x  =  20s 


T.     =  691. 8  K 
in,o 
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the  feedback  control  given  by  Eqs.  (2-16)  through 

(2-18).   The  control  problem  is  how  the  inlet  temperature 

should  be  manipulated  whenever  the  outlet  conversion 

decreases  to  a  certain  level  from  the  desired  conversion 

(bandwidth) .   Note  that  the  intrinsic  rate  of  deactivation 

rp  is  used  here  only  to  generate  the  process  response. 

In  the  feedback  control,  r   is  treated  as  an  unknown  and 

P 

only  the  on-line  measurements  (simulated  temperature  and 
concentration)  are  used  to  manipulate  the  inlet  temperature. 
The  constraints  are  the  maximum  reactor  temperature  allowed 
and  the  final  time  at  which  the  catalyst  is  regenerated  as 
the  time  when  the  final  bandwidth  becomes  one-tenth  of  the 
initial  bandwidth.   These  constraints  are  also  given  in  Table 
3-2  along  with  the  desired  outlet  concentration  (conversion) 
and  the  reactor  size.   The  initial  inlet  temperature  for 
fresh  catalyst  is  the  one  corresponding  to  the  specified 
reactor  size  (t)  and  outlet  conversion,  which  is  691. 8K. 
Since  the  film  heat  transfer  coefficient  is  quite  small  for 
the  example  problem  and  thus  the  value  of  q  is  relatively 
large,  the  proportionality  factor  6  in  Eq.  (2-16)  was  set  at 
0.95. 

The  behavior  of  the  model  reactor  resulting  from  the 
feedback  control  is  shown  in  Figure  3-1.   Initially,  the 
inlet  temperature  is  at  its  initial  value  and  the  outlet 
temperature,  which  is  the  maximum  reactor  temperature  in 
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this  case,  is  well  below  the  maximum  temperature  allowed. 

Due  to  the  deactivation,  however,  both  the  conversion  and 

the  outlet  temperature  decrease  with  time  and  the  conversion 

eventually  reaches  the  lower  bound  given  by  the  bandwidth, 

triggering  an  adjustment  of  the  inlet  temperature.   The  new 

inlet  temperature  (T.),  calculated  by  the  algorithm  of 

Eqs.  (2-16)  through  (2-18)  is  seen  to  bring  the  conversion 

back  to  the  desired  level,  the  actual  value  being  6  9.96%  as 

opposed  to  the  desired  value  of  70%.   A  bisection  method 

was  used  for  the  calculation  of  (T.  ) ..   When  the  outlet 

in  j 

temperature  reaches  the  allowed  maximum,  which  occurs  at 

4 
around  t  =  2x10  s,  the  conversion  can  no  longer  be  brought 

back  to  the  desired  level  due  to  the  temperature  constraint. 
Consequently  the  bandwidth  decreases  as  time  increases  and 
eventually  the  bandwidth  becomes  less  than  one-tenth  of 
the  original  bandwidth,  resulting  in  the  reactor  shut-down 
for  catalyst  regeneration  according  to  the  constraint  im- 
posed.  In  practice,  however,  the  conversion  may  be  allowed 
to  decrease  below  60%  when  the  maximum  temperature  is  reached. 
In  the  example,  the  average  conversion  is  essentially  main- 
tained at  the  middle  of  the  bandwidth,  which  is  65%. 
§3-2  Experimental  Verification 

We  have  shown  the  feasibility  of  control  policies  for 
the  model  reaction  system.   The  theory  has  also  been  put 
to  test  using  a  shell-and-tube  reactor.   Experimental  results 
obtained  and  comparisons  with  the  theory  are  presented  in  this 
section. 
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§3-2-1  Catalyst  Preparation 

The  catalyst  used  in  this  study  of  ethylene  oxidation 
was  silver  supported  on  fused  alumina.   The  quarter-inch- 
diameter  alumina  obtained  from  Norton  Chemical  Company 

(Catalog  No.  SA-5202)  was  crushed  into  8-9  mesh  (0.20- 

0.24cm)  with  a  rod  miller,  washed  with  distilled  water, 

and  then  dried  in  the  furnace  at  37  3 K overnight .   These 

particles  were  porous  with  0.7-1.3  m  /g  BET  surface  area 

3 
and  1.9  g/cm  bulk  density.   The  alumina  was  impregnated 

with  an  aqueous  solution  of  lactic  acid  and  silver  oxide 
at  368  K  for  90  minutes.   The  solution  contained  silver 
oxide,  85%  lactic  acid  and  distilled  water  in  the  ratio  of 
1:2:2  by  weight.   After  impregnation  the  excess  liquid 
was  removed  and  the  particles  were  heated  in  a  furnace  for 
16  hours  at  643  K  to  decompose  the  silver  lactate  to  metal- 
lic silver.   The  catalyst  was  further  stabilized  by  oxy- 
genating at  573  K  for  100  hours.   Since  the  color  of  these 
pellets  was  not  uniform,  only  those  with  lighter  color 
were  picked  for  the  experiment.   The  weight  gain  indicated  that 
the  catalyst  contained  12.4%  silver.   The  oxygen  chemisorp- 
tion  at  423  K  was  36  \ii/g. 
§3-2-2  Experimental  Apparatus 

The  reactor  used  for  ethylene  oxidation  reactions 
was  shell-and-tube  type  as  shown  in  Fig.  3-2.   The  shell- 
and-tube  section  was  6  0.1  cm  long  with  outside  diameter 
(o.d.)  of  3.81  cm  for  the  shell  and  1.19  cm  for  the  tube. 
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One  thermal  well  (A,  0.794  cm  o.d.)  ran  axially  into  the  tube, 
another  (B,  0.0318  cm  o.d.)  into  the  shell.   Silver  catalyst 
was  packed  into  the  annular  section  between  the  tube  and 
thermal  well  A.   The  shell-and-tube  consisted  of  three 
segments  (from  left  to  right  in  Fig.  3-2):   preheating  zone 
(20.3  cm),  reaction  zone  (35.6  cm)  and  outlet  zone  (5.1  cm). 
Glycerin  was  used  as  coolant  which  circulated  through  the 
shell  side  using  a  Haake  HT22  temperature  controller/circu- 
lator.  A  total  of  twenty  temperatures  were  measured  as 
shown  in  Fig.  3-2:   eighteen  in  thermal  well  A  and  two  in 
B. 

All  temperature  points  (thermocouples)  were  connected 
to  a  rotary  selector  switch  which  was  in  turn  connected  to 
a  temperature  indicator.   The  feed  gases  were  introduced  into 
the  reactor  at  the  desired  mass  flow  rates  and  composition 
controlled  by  a  multiple  channel  electronic  mass  flow  con- 
troller.  Two  pressure  gauges  were  installed,  one  at  the 
inlet,  the  other  at  the  outlet.   The  outlet  composition  was 
analyzed  by  a  Tracor  550  gas  chromatograph.   The  coolant 
flow  rate  was  set  high  enough  to  maintain  an  essentially 
isothermal  reaction  zone.   The  reactor  pressure  was  held 
constant  at  1  atm. 
§3-2-3   Determination  of  Kinetics 

Ethylene  oxidation  involves  a  triangle  reaction  net- 
work (Spath  and  Handel,  1974;   Dettwiler  et  al. ,  1979) 
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1 


C2H4 
(A) 


*r    C2H4° 

/     (B) 


C02  +  H20 
(C) 

The  oxidation  of  ethylene  oxide  (reaction  path  3)  is 

significant  only  at  temperatures  above  473  K  (Dettwiler 

et  al.,  1979;  Akella,  1983).   Since  the  experiments  were 

conducted  in  the  temperature  range  between  448  and  4  68  K, 

only  epoxidation  and  combustion  of  ethylene  (reaction  paths 

1  and  2)  need  to  be  considered.   A  material  balance  in  an 

isothermal  fixed-bed  can  be  written  as 


dz 


=  -  h1(r1  +  3r2)W 


=  F 


z=0 


Ao 


(3-1) 


dF 


B 


dz 


=  h^W 


B 


z=0 


=  FBQ  =  0  (3-2) 


where  FA  and  FB  are  the  molar  flow  rates  for  ethylene  and 
ethylene  oxide,  respectively,  and  w  is  the  weight  of  the 
catalyst.   Here  r^  and  r2  are  the  intrinsic  rates  for  the 
reaction  paths  1  and  2.   The  mechanism  of  ethylene  oxidation 
is  still  controversial.   Under  the  condition  of  constant 
oxygen  partial  pressure,  however,  the  following  Langmuiur- 
Hinshelwood  rate  expressions  are  reported  by  most  investiga- 
tors (Dettwiler  et  al.,  1979): 


30 


k,k  P, 
r,  =   ,  !„**  (3-3) 


1     l+k  p^ 
q^A 


k  k  p 
ro  =  uA —    (mol/s-g-cat.)  (3-4) 

2    i+KqPA 

where  p  is  the  partial  pressure  of  ethylene  and  k   is 
an  adsorption  equilibrium  constant.   Since  the  catalyst 
particles  are  small,  there  is  no  diffusion  limitation  and 
the  intrinsic  rates  are  essentially  the  same  as  the  glo- 
bal rates. 

Let  x  and  y  be  the  conversion  and  yield  in  the  reactor, 

i'e-  X  =  (FAo  "  FA)/FAo  and  y  =  (FB  "  FBo)/FAo-   Then' 
Eq.  (3-1)  can  be  rendered  dimensionless  with  the  aid  of 

Eqs.  (3-3)  and  (3-4): 


dx    _W_  WBk2)  kqPyko(l-x) 
dz  "  Vko      1+VyAo(1-x) 


x 


=  0 
z=0 


(3-5) 


where  FQ  is  the  molar  flow  rate  of  reactant  mixture  and 
y^Q  is  the  mole  fraction  of  ethylene,  both  at  the  inlet  of 
the  reactor.   The  ideal  gas  law  has  been  applied  to  obtain 
Eq.  (3-5).   The  reactor  activity  factor  for  the.reaction  path  1, 
H1  can  be  obtained  by  rearranging  Eq.  (3-5)  and  integrating 
from  the  reactor  inlet  to  the  outlet: 


Hl  " 


r1       kaPykox-£n(1_x) 

hidz  =    (5.4?  )k  oW/F  (3-6> 

1    2.      q^    o 


0  l    ^  q 

where  X  is  the  conversion  at  the  outlet 
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Dividing  Eq.  (3-1)  by  (3-2)  and  nondimensionalization  gives 


dx  =  x  +Bf2  =  1   +B^2 
dy        r,        k 


=  0 


(3-7) 


y=o 


Since  the  right  hand  side  of  Eq.  (3-7)  is  constant  through- 
out the  isothermal  reactor,  Eq.  (3-7)  can  be  integrated  from 
the  reactor  inlet  to  the  outlet  to  give 


*  =  1  + 
Y 


(3-8) 


where  Y  is  the  yield  at  the  outlet. 

Equations  (3-6)  and  (3-8),  which  are  the  integrated 
versionsof  Eqs.  (2-26)  and  (2-27) ,  can  be  used  for  both 
the  feedback  control  described  in  Section  §2-3-3  and  the 
determination  of  kinetics.   Before  proceeding  to  the  control 
policies,  consider  the  kinetics  first.   It  is  convenient 
to  choose  the  activity  of  the  catalyst  at  the  time  kinetic 
data  are  taken  as  the  reference  state  activity.   Then, 
both  h,  and  6  are  unity  by  definition,  and  Eqs.  (3-6) 


and  (3-8)  can  be  rearranged  as 


and 


[£n  l-x] 

-k  p 

q 

y»  x 

J  Ao 

+ 

KiVw 

Fo 

X1 

i  Y 

k2    =   kl 

Y 

-  i 

(3-9) 


(3-10) 


It  has  been  shown  (Akella,  1983)  that  the  selec- 
tivity, defined  as  S  =  Y/X,  is  fairly  independent  of  ethylene 
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partial  pressure  but  is  dependent  on  the  reaction  tem- 
perature.  This  can  also  be  seen  from  Eq.  (3-8)  where  the 
left  hand  side  is  1/S  and  the  right  hand  side  depends  only 
on  temperature.   Thus,  at  any  given  temperature,  X/Y  is 

constant,  and  Eqs .  (3-9)  can  be  used  to  obtain  k   and  k  k 

q      1  q 

by  a  linear  regression  of  ( £n  ^r^)  on  WA  x> •   The  rate 
constant  k2  can  then  be  calculated  from  Eq.  (3-10).   By  de- 
finition and  the  stoichiometric  relationships,  the  conversion 
and  yield  at  the  reactor  outlet  can  be  related  to  the  inlet 
and  outlet  gas  compositions  as 

X   =  1      ~  yT-  (3-11) 

JAo 
y" 

Y  -  X  "  #  (3-12> 

yAo 

where  y^Q  and  y^  are  the  mole  fractions  of  ethylene  at  the 
inlet  and  outlet,  respectively,  y«  is  the  mole  fraction  of 
carbon  dioxide  at  the  outlet. 

We  have  shown  that  at  a  given  temperature,  the  rate  con- 
stants, k   and  k„,  and  the  absorption  equilibrium  constant,  k  , 

^  q 

can  be  determined  from  the  inlet  and  outlet  gas  compositions, 
YAo'  yAl  and  yCl'  usin9  Eqs-  (3-9)  through  (3-12).   The  kinetic 
data  for  ethylene  oxidation  reactions  are  given  in  Table  3-3. 
The  reactor  was  maintained  isothermal  by  controlling  the  coolant 
flow  rate  at  a  high  level.   Due  to  the  low  ethylene  concentra- 
tion and  the  high  coolant  flow  rate,  the  heat  generated  by  the 
reactions  was  removed  immediately  by  the  coolant.   The  oxygen  com- 
position in  the  feed  gas  was  maintain  at  20%  for  all  runs,  with 
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helium  as  the  diluent.   Given  in  Table  3-4  are  the  results 

of  the  linear  regression  of  (Jin  ^~)  on  (y^QX)  for  each 

temperature.   The  values  of  (in  -^)    calculated  from  the 

regression  model,  (in   ^-)    ,   are  included  in  the  table 

A  cal 
for  comparison.   The  percent  error  in  the  table 

was  calculated  from 

1 


Error  = 


In 


!"X   -  -  1 


(£n  T^> 

1  X  Cal 


x  100%  (3-13) 


The  constants  k±   and  kg  were  calculated  from  the  slope  and 

intercept  of  the  regression  line,  while  k2  was  calculated 

from  Eq.  (3-10) .   The  values  of  these  rate  constants  are 

given  in  Table  3-4  along  with  the  coefficients  of  correlation. 

Note  that  the  average  selectivity,  £,  given  in  the  last  column 

of  Table  3-3  was  used  to   calculate    kn  and  k 

1     q" 
It  was  assumed  that  the  temperature  dependence  of  the 

rate  and  equilibrium  constants  are  all  of  the  Arrhennius 
form: 

k.  =  A.  exp(-  E./RgT)    i=l,  2,  q  (3-14) 

The  activation  energies  and  pre-exponential  factors  were 
obtained  from  linear  regressions  of  In  k.  on  1/T.   The  re- 
sults are  given  in  Table  3-5. 

The  experimental  results  show  that  the  kinetics  of 
ethylene  oxidations  reactions  can  be  well  represented  in 
the  form  of  Eqs .  (3-3)  and  (3-4)  with  the  following  rate 
constants  and  equilibrium  constant: 
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k±   =  3.618  x  10~2  exp(-7098/T)  (mol/s . g-cat ) 

k2  =  3.038  x  10_1  exp(-8279/T)  (mol/s . g-cat) 

kg  =  1.942  exp  (1641/T)  (1/atm)  (3-15) 

0.04  atm   <  P   <  0.08  atm 
Ao 

448  K  <_   T  <    468  K 

The  comparisons  made  in  Table  3-6  between  the  experi- 
mental conversions  and  those  predicted  by  the  kinetic 
model  show  very  good  agreements,  the  maximum  error  being 
1.9%.   The  same  can  be  concluded  for  the  yield  and  selec- 
tivity.  It  is  seen  from  the  table  that  the  maximum  errors 
are  2.3  and  3.4%  respectively  for  the  yield  and  selectivity. 

The  experiments  for  the  kinetic  data  lasted  9  hours. 
An  additional  run  carried  out  right  after  the  period  with 
the  same  conditions  as  in  Run  No.  2  gave  almost  the  same 
conversion  and  yield  as  shown  in  Fig.  3-3  at  t  =  o  and  t  =  9 
hr.   This  is  a  good  indication  that  no  significant  deacti- 
vation had  occurred  during  the  kinetic  experiments. 
§3-2-4  Experimental  Verification  of  the  Control  Policy 

The  shell-and-tube  reactor  at  the  end  of  the  kinetic 
experiments  was  allowed  to  run  continuously  to  verify  the 
feedback  control  policy  presented  in  Chapter  2.   The  feed 
rate  was  controlled  at  100  ml/min  with  the  ethylene, 
oxygen  and  helium  concentrations  of  6 ,  20  and  74%,  respectively, 
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The  temperature  was  kept  constant  at  44  8  K.   The 
outlet  gas  compositions  were  analyzed  daily  and  the  corres- 
ponding conversion  and  yield  were  calculated.   The  control 
policy  was  that  of  maintaining  the  yield  of  ethylene  oxide 
at  8%  in  a  piecewise  manner.   The  bandwidth  was  set  in  such 
a  way  that  an  adjustment  was  made  whenever  the  yield  de- 
creased to  6.5%.   The  temperature  adjustments  and  the 
corresponding  reactor  behavior  over  a  period  of  35  days 
are  shown  in  Fig.  3-3.   The  first  temperature  adjustment 
occurred  at  t  =  120  hr.  at  which  time  the  yield  dropped 
down  to  6.5%.   This  triggered  an  adjustment  of  the  tem- 
perature in  accordance  with  the  specification  of  the  band- 
width . 

Consider  the  temperature  adjustment,  in  particular 

the  procedures  of  calculating  the  new  temperature  which 

will  bring  the  yield  back  to  the  desired  level.   The  value 

of  6  required  to  calculate  H,  from  Eq.  (3-6)  is  obtained 

directly  from  the  experimental  values  of  the  conversion  and 

yield  at  the  current  temperature.   This  in  turn  is  used  in 

Eq.  (3-6)  to  calculate  the  current  reactor  activity 

factor  for  the  reaction  path  1,  Hlc.   Rewriting  Eq.  (3-6) 

for  the  new  temperature  and  using  the  relationship  H  =  H 

c      lc    In 

yields 

r  kgPyAo  xd  -  ^(1  -  Xd) 


lc 


(k   +  6k_)  k  PW/F 
1      2    q     o 


*  Jn 


(3-16) 
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Here,  the  subscript  n  denotes  the  new  temperature  and  X 

d 

is  the  desired  conversion.   Since  local  activities  change 
little  during  the  short  period  of  temperature  adjustment, 
B,  by  definition,  also  changes  little.   This  means  that 
the  same  3  obtained  from  Eq.  (3-8)  can  be  used  in  Eq. 
(3-16)  to  calculate  the  new  temperature.   If  the  conversion 
Xd  is  the  one  to  be  controlled,  then  Eq.  (3-18)  can  be 
used  to  calculate  Tn#  since  it  is  the  only  unknown  in  Eq. 
(3-18).   In  this  case,  however,  the  yield  Y   is  to  be  con- 
trolled.  Hence,  we  need  an  additional  relationship  between 
Xd  and  Yd*   This  relationship  can  be  obtained  by  rewriting 
Eq.  (3-8)  for  the  new  temperature 


^=1  +  6 
d 


k 


2 


kl 


(3-17) 

n 


This  relationship  was  used  for  the  calculation  of  the  new 
temperature  in  the  experiment. 

The  experiments  were  conducted  over  a  period  of  8  84 
hours.   As  shown  in  Fig.  3-3,  the  first  adjustment  was 
made  at  t  =  120  hr.   The  temperature  was  raised  to  454  K, 
which  brought  the  yield  to  7.95%  for  the  desired  yield  of  8%. 
Two  additional  temperature  adjustments  were  made  at  t  =  336 
hr.  and  t  =  6  24  hr.   The  experimental  values  of  the  yield  of 
ethylene  oxide  resulting  from  the  temperature  adjustements  were 
both  7.92%.   The  behavior  of  conversion  and  selectivity  are 
also  shown  in  Fig.  3-3.   It  is  interesting  to  observe  that 
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the  selectivity  increases  at  448  K,  remains  almost  constant 
at  454  K  and  then  decreases  at  460  and  466  K.   This  behavior 
can  be  explained  in  terms  of  H1  and  H  .   As  shown  in  Fig. 
3-4,  H1  and  H2  decrease  with  different  speeds,  resulting  in 
the  observed  selectivity  behavior.   The  results  given  in 
the  figure  show  excellent  agreements  between  the  theory  and 
the  experiment. 

The  experimentally  determined  values  of  H.  and  H  were 
fitted  to  an  nth  order  deactivation  kinetics  using  a  non- 
linear regression.   The  results  are 

dH 

dt~  =  "  kD   H^1 =-5.018  x  105  exp  (-8879/T) hJ " 87 

dH2 

dt~  =  "  kD   H22  =  "  1'936  x  10~8  exP  (5707/T)H^19 
2  2 

The  experimental  values  are  plotted  in  Fig.  3-4  along  with 
the  calculated  values,  which  are  shown  as  solid  curves. 
§3-3  Concluding  Remarks 

We  have  shown  in  Section  §3-1  the  feasibility  of  the 
feedback  control  policies  by  simulating  an  adiabatic  fixed- 
bed  where  a  single  reaction  affected  by  diffusion  and  poi- 
soning takes  place.   The  feedback  control  policies  have  also 
been  put  to  test  with  a  laboratory  isothermal  fixed-bed 
where  multiple  reactions  affected  by  sintering  take   place. 
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The  results  in  Fig.  3-1  and  3-3  show  that  one  can  indeed 
maintain  the  conversion  or  yield  at  the  desired  level  with 
good  accuracy.   It  is  noted  in  this  regard  that  the  con- 
trol policies  can  also  be  used  to  control  the  conversion 

or  yield  in  any  desired  manner,  since  X  or  Y,  can  be 

d     a 

changed  each  time  the  temperature  is  raised. 

The  accuracy  of  the  feedback  control  depends  mostly  on 
the  on-line  measurements  (X  and  Y)  but  slightly  on  the  reference- 
state  kinetics  determined.   For  example,  at  t  =  624  hr., 
(as  shown  in  Fig.  3-3  when  T  is  adjusted  from  460  to  466  K) , 
5%  error  in  the  determination  of  yield  results  in  1.7K  dif- 
ference in  the  new  temperature  calculated,  while  5%  error  in 
activation  energy  for  reaction  path  1  (E^  leads  to  only  0.3  K 
difference.   In  calculating  new  temperature  by  use  of  Eq.  (3-16), 
the  on-line  measurements  of  conversion  and  yield  are  used  only 
on  the  left  hand  side  of  the  equation,  while  the  activation 
energy  E1  is  used  on  both  sides.   Thus,  any  error  in  E   is 
likely  to  be  cancelled  out  from  both  sides  of  the  equation. 
This  is  a  good  indication  that  the  feedback  control  is  appli- 
cable even  when  some  uncertainty  exists  in  the  kinetic  model. 


CHAPTER  4 
OPTIMAL  PIECEWISE  CONTROL 

The  control  policies  discussed  in  the  previous 
chapters  are  all  piecewise  type,  i.e.  the  control  actions 
are  taken  at  discrete  intervals.   The  optimal  control 
problems  for  reactors  affected  by  catalyst  deactivation  have 
been  the  subject  of  several  studies  (Chou  et  al.,  1967; 
Szepe  and  Levenspiel,  1968;  Crowe,  1970),  and  a  constant 
conversion  policy  has  been  proved  to  be  optimal  for  a 
single  irreversible  reaction  with  concentration-indepen- 
dent deactivation.   These  studies  deal  only  with  the  con- 
tinuous control,  which  is  essentially  the  piecewise  control 
with  an  infinite  number  of  steps.   However,  the  feasibility 
of  extending  the  optimal  continuous  control  policies  to 
piecewise  control  schemes  with  a  finite  number  of  steps 
remains  questionable.   In  this  chapter,  a  solution 
method  for  the  optimal  piecewise  control  is  developed  and 
a  comparison  between  continuous  and  piecewise  control  is 
made. 
§4-1   Problem  Statement  and  Solution  Method 

Consider  an  isothermal  fixed-bed  where  a  single  irre- 
versible reaction  affected  by  catalyst  deactivation  takes 
place.  The  temperature  is  to  be  adjusted  in  steps  from  T,  to 
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Tn  as  shown  in  Fig.  4-1.   The  pseudo  steady-state  balance 
equations  can  be  written  as 


dt  =  g(Tlf  T2,. . .T^,h)  ;  h 


t-o  =  x    (4-1} 

X  =  f(T1#  T2,...Tn  ,  h)  (4-2) 

where  h  is  the  activity  of  the  catalyst,  X  is  the  conver- 
sion at  the  outlet  of  the  reactor,  and  g  represents  the 
dependence  of  deactivation  rate  on  the  temperatures  (T) 
and  activity  (h) .   Equation  (4-2)  is  obtained  by  integrating 
the  mass  balance  equation  from  the  inlet  to  the  outlet  of  the 
reactor. 

The  optimization  problem  is  to  maximize  the  total 

amount  of  feed  converted  over  a  fixed  total  reaction  time   t 

n' 

by  choosing  the  best  set  of  temperatures,  [T  ,  T      t  1 

1    2' "  " " '  n 


That  is 


J  =  max  J       i  =  1,  2,...n  (4-3) 

T. 

l 


with 

t  t 

1    f  n       1        n      (   i 

Xdt  (4-4) 


n 


t 

t . 

n 

1 

n 

C   i 

Xdt 

-■£, 

0 

t 

i=l 

n 

t.    , 

i-I 

and  subject  to 


T  .  <.  t.    T 
mm     i  —  max 


t 
At  =  ti   ~   t-L-i  =  constant  =  -2  i  =  lf    2, 


n 

(4-5) 
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Equation  (4-1)  can  be  solved,  analytically  or  nu- 
merically, to  give 

h  =  h(T1,  T2,...Tn,  t)  (4-6) 

Integrating  Eq.  (4-4)  with  the  aid  of  Eqs.  (4-2)  and  (4-6) 
yields 

J  =  J(T. ,  T_,...T  )  (4-7) 

1   2.  n 

Thus,  the  optimal  piecewise  temperature  control 

problem  of  Eqs.  (4-3)  through  (4-5)  is  a  multivariable 

optimization  problem  searching  for  n  variables  (T. ,  T„  ,  .  .  . 

T  ) ,  rather  than  a  variational  calculus  problem  searching 
n 

for  a  continuous  function  (T(t)).  The  latter  is  often  en- 
countered in  problems  dealing  with  optimal  continuous  con- 
trol. The  stationary  point  of  J  can  be  obtained  either  by 
applying  the  necessary  conditions 

|f-  =  0       i  =  1,  2,...n         (4-8) 
i 

and  then  solving  these  n  nonlinear  equations  numerically, 
or  by  direct  numerical  search  (Beveridge  and  Schechter, 
1970) .   It  is  worth  noting  that  even  for  the  simplest  first 
order  reaction  and  deactivation,  the  expression  for  J  is 
still  very  complex  and  the  partial  derivatives  in  Eq. 
(4-8)  are  difficult  to  obtain  analytically.   Therefore, 
direct  numerical  search  is  a  better  method  for  obtaining 
the  optimal  temperature  policy,  T, .  T„...,T  .  The 
optimal  temperature  policy  leads  directly  tothe  optimal  conver- 
sion policy  through  Eq.  (4-2)  with  the  aid  of . Eq .  (4-1). 
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§4-2   Optimal  Piecewise  Control  Policies 

As  an  illustration,  the  optimal  piecewise  control 
problem  of  Eqs.  (4-1)  through  (4-5)  was  solved  for  a  model 
reaction  system  of  first  order  irreversible  reaction  with 
first  order  deactivation  (Table  4-1) .   The  multivariable 
optimization  problem  was  solved  by  direct  numerical  search 
using  the  ZXMWD  subroutine  from  International  Mathematical 
and  Statistical  Libraries  (IMSL) .   The  ZXMWD  subroutine 
is  used  for  obtaining  the  global  minimum  or  maximum  for  a 
multivariable  function  with  constraints.   The  constrained 

maximization  problem  of 

* 

J   =    max    J(T  ,  T  , ...T  )  (4-9) 

Tm        x     ^       n 

1     n 

with 

T.   <  T.<   T      i  =  1,  2,  . . .  n 

min  -   i  -   max        '   ' 

is  transformed  in  ZXMWD  subroutine  to  an  equivalent,  but 
unconstrained  problem 

*  ?- 

J      =      max  J[T    .    +(T         -   T    .     )Sin   T,,...,T    .       + 

=  =  mm        max        mm  1'         '   mm 

1  n 

(Tmax  "  Tmin)Sin2  n]  (4"10) 

where  T^  is  the  transformed  function  of  T^and  is  now  uncon- 
strained.  With  this  transformation,  each  possible  global 
maximum,  including  any  on  the  boundary,  is  transformed  into  a 
local  maximum  (Box,  196  6) . 

Shown  in  Figure  4-2  are  the  optimal  piecewise  control 
policy  for  the  model  reaction.   The  optimal  temperature 
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Table  4-1 Equations  and  Parameters  for  a  Model  Reaction 

System  for  Optimal  Piecewise  Control 


dh 
dt 


-v 


hlt=0=  1 


X  =  1  -  expf-xk-h) 


kD  =  AD  exP 


(      Gr 


k,  =  A1  exp 


f      Gn 


t  t 

T(t)  =  T     for   (i-l)-S  <L  t  <  ±-2  ;   i  =  i  2     n 


T  .   <.  T.   :  T 
mm    i  ~  max 


i  =  1,2  , . . .n 


AD  =  5.00  x  106  hr-1 
GD  =  1.25  x  104  K 
tn  =  2000  hr 


Ax  =  5.00  x  106   S_1 


G1   =   1.00  x  10 
T  =  20  S 


K 


mm 


max 


n 


=  473  K 

=  573  K  in  Figs.  (4-2),  (4-3)  and  (4-5) 

=  673  K  in  Figs.  (4-4),  (4-6)  and  (4-7) 

=  10     in  Figs.  (4-2)  through  (4-6) 

=  varies  in  Fig.  (4-7) 
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policy  is  to  step  the  temperature  up  from  539.9  K  until  it 
reaches  Tmax  (573  K)  at  t  =  1400  hr.  and  then  to  keep  it  at 
that  level  as  shown  in  the  figure.   The  corresponding  con- 
version zigzags  with  increasing  bandwidth  until  the  tempera- 
ture reaches  its  upper  limit.   Due  to  catalyst  deactivation 
the  conversion  decreases  in  each  time  interval  while  the 
temperature  remains  constant.   It  is  of  interest  to  note 
that  each  temperature  increase  gives  a  conversion  peak 
which  is  slightly  higher  than  the  previous  one,  i.e. 

Xn(ti)  >    Xn(ti  1>  for  T  <  T 
n  ±  n   1--L  max 

where  n  denotes  the  new  conversion  immediately  following  the 
temperature  adjustment.   In  other  words,  the  optimal  piece- 
wise  policy  is  an  increasing  conversion  policy,  rather  than  a 
constant  conversion  policy  which  is  the  optimal  for  continuous 
control  (Crowe,  1970).   The  optimal  continuous  control  policy 
for  the  same  model  reaction  is  shown  in  Fig.  4-3  for  compari- 
son.  It  is  obtained  by  one-dimensional  search  of  the  optimal 
initial  temperature  with  the  aid  of  the  constant  conversion 
policy.   it  can  be  seen  from  the  figure  that  the  optimal  con- 
tinuous policy  is  a  constant  conversion  policy  followed  by  a 
constant  temperature  policy  at  Tmax.   Comparison  between  Fig. 
4-2  and  4-3  shows  that  the  temperature  in  the  continuous  con- 
trol starts  1.3  K  lower  than  in  the  piecewise  scheme  and  reaches 
the  upper  limit  42  hours  later.   It  can  also  be  seen  that  the 
conversion  in  the  continuous  control  is  lower  than  the  original 
conversion  in  the  piecewise  control,  but  stays  in  the  zigzag 
zone  of  the  piecewise  scheme.   The  corresponding  activity  be- 
havior is  also  plotted  in  Figs.  4-2  and  4-3. 
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Another  major  difference  between  the  piecewise  and  the 
continuous  control  is  that  although  the  temperature  in  the 
optimal  continuous  control  always  ends  at  its  upper  limit 
for  the  model  reaction  studied  here  (Crowe,  1970)  ,  the  tem- 
perature in  the  optimal  piecewise  control  may  end  at  a  tem- 
perature below  the  upper  limit.   This  is  shown  in  Fig.  4-4 
where  the  final  temperature  is  621.3  K,  far  below  the  upper 
limit,  673  K.   However,  the  final  temperature  will  approach 
673  K  as  the  number  of  control  steps,  n,  increases. 

A  comparison  of  performance  indices  calculated  from 
Eq.  (4-4)  for  seven  different  control  policies  is  given  in 
Table  4-2  with  Tmax  =  573  K  and  673  K.   The  performance  in- 
dices calculated  with  the  optimal  piecewise  temperature 

* 

control  policy,  T,  are  listed  in  Entry  1.   In  Entries  2  and  3, 

the  temperature  control  policies  are  1  K  higher  and  lower, 

* 

respectively,  than  T  ,  but  still  within  the  interval 

[Tmin'  Tmax] -   The  contro1  policy  Tc  in  Entry  4  is  obtained 
by  solving  the  same  optimization  problem  with  an  additional 
constraint  of  constant  conversion.   Note  that  by  applying 
the  constant-conversion  constraint,  the  multivariable  optimi- 
zation problem  of  Eq.  (4-9)  is  simplified  to  a  single  vari- 
able optimization  problem 

Jc  =  max  JC(T1)  (4-11) 

Tl 
since  T2,  T3,...  and  Tn  are  fixed  when  T   is  specified. 

Here  subscript  c  denotes  the  performance  index  obtained 

under  the  constant-conversion  constraint.   The  performance 
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indices  obtained  from  Eq.  (4-11)  are  listed  in  Entry  4.   In 

Entries  5  through  7 ,  the  temperature  policies  are  all 

constant  with  respect  to  time.   It  can  be  seen  from  the  table 

* 
that  T   (the  optimal  temperature  policy  shown  in  Fig.  4-2) 

is  indeed  the  optimal  temperature  control  policy,  since  all 
the  other  policies  result  in  smaller  values  of  the  performance 
index. 

The  single-variable  optimization  problem  of  Eq.  (4-11) 
was  solved  numerically  using  Fibonacci  search  (Beveridge 
and  Schechter,  1970).   The  results  are  shown  in  Fig.  4-5 
f°r  Tmax  =  573Kand  in  Fig.  4-6  for  Tmax  =  673  K.   It  can  be 
seen  from  the  figures  that  the  conversion  is  brought  back 
to  the  initial  conversion  each  time  the  temperature  is  ad- 
justed. 

§4-3   Concluding  Remarks 

A  comparison  between  the  optimal  performance  indices 
obtained  with  and  without  constant-conversion  constraint  is 
shown  in  Fig.  4-7  as  a  function  of  the  number  of  control  ac- 
tions (n) .   it  is  seen  that  as  n  approaches  infinity,  J*  and 

* 
Jc  approach  the  same  asymptotic  value  (0.56556  for  the  model 

reaction  with  T^  =  673  K)  .   This  indicates  that  the  con- 
tinuous control  (which  corresponds  to  an  infinite  number  of 
piecewise  control  actions)  is  superior  to  the  piecewise 
control  in  terms  of  the  performance  index.   On  the  other 
hand,  J  >  99.5%  when  n>21  as  shown  in  Fig.  4-7,  implying 
that  the  piecewise  control  is  almost  as  good  as  the  con- 
tinuous control  if  the  number  of  control  actions  is  large 
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enough.   In  practice,  the  temperature  adjustments  are  made 
intermittently  because  of  the  slow  nature  of  catalyst  de- 
activation in  fixed-beds.   From  the  fact  that  the  curve 
for  J  lies  above  the  curve  for  J  ,  it  can  be  concluded 
that  the  optimal  piecewise  control  policy  is  an  increasing 
conversion  policy  as  shown  in  Fig.  4-2  and  4-4.   However, 
the  constant-conversion  policy  can  be  used  when  number  of 
control  steps  is  large,  say,  30. 

The  problem  formulation  and  solution  method  presented 
in  Section  4-1  is  neither  restricted  to  an  isothermal 
reactor  nor  to  a  single  irreversible  reaction.   For  example, 
one  can  consider  a  reversible  reaction  taking  place  in  an 
adiabatic  fixed-bed.   In  this  case,  the  reactor  tempera- 
tures, T^s,  in  Eqs.  (4-1)  through  (4-5)  should  be  replaced 

by  the  temperatures  at  the  inlet  of  the  reactor,  (T.  ) . 's 

in  i 

Equation  (4-2)  might  then  become  implicit  in  X  because  of 
the  reversible  reaction  and  the  nonisothermal  reactor. 
However,  X  can  still  be  solved  numerically.   Thus,  the 
problem  remains  a  multivariable  optimization  problem  which 
can  be  solved  by  application  of  an  appropriate  numerical 
method. 


CHAPTER  5 
OPTIMAL  CONTINUOUS  CONTROL  FOR  COMPLEX  REACTION 

In  the  experiment  described  in  Section  §3-2,  we  con- 
trolled the  yield  of  ethylene  oxide  at  a  constant  level  in 
a  piecewise  manner,  and  verified  the  feasibility  of  the 
feedback  control  policy  presented  in  Chapter  2. 

However,  the  constant  yield  policy  used  is  not  an 
optimal  policy.   It  is  the  purpose  of  this  chapter  to  ob- 
tain an  optimal  control  policy  for  complex  reactions  which 
are  parallel  with  the  Langmiur-Hinshelwood  rate  expressions. 
We  consider  here  only  the  optimal  continuous  control  policy 
since  we  have  shown  in  the  previous  chapter  that  one  can 
use  the  optimal  continuous  control  policy  for  piecewise  con- 
trol as  long  as  the  number  of  control  steps  is  large  enough. 
§5-1   Problem  Statement 

Consider  a  parallel  reaction  network  with  the  Langmiur- 
Hinshelwood  rate  expressions  for  an  isothermal  fixed-bed  such 
as  the  ethylene  oxidation  reaction  described  in  Section  §3-2-3 
Assume  the  same  active  sites  are  responsible  for  both  reac- 
tions, i.e.,  h1  =  h2  =  h  and  6  =  1.   Equation  (3-5)  and  (3-7) 
can  be  rewritten  as 

dx   Th(Kl  +  V  (1_x) 


dz       1  +  K  (1-x) 


z=o=0  <5-D 


and 


dy       kx 


y=o  =  °  (5-2) 


62 


63 


where 


x  =  WP/F  ,  K,  =  k,k  ,  K  =  k0k  ,  and  K  =  P^  k 
o    1     lq    2    2q'       q    Aoq 


The  rate  of  decay  of  the  catalyst  activity  can  be  expressed 
as 


dh      T,  ,  m  . 

dt   =  "  KDh         '   h 


z=o=  1  <5"3> 


where  K  is  the  deactivation  rate  constant.  Integrating 
Eqs.  (5-1)  and  (5-2)  from  the  inlet  to  the  outlet  of  the 
reactor  yields 

K  X  -  £n(l-X)  =  Thfl^  +  K2)  (5-4) 

Kl 
Y  =   Kpn^  X  (5-5) 

where  X  and  Y  are  the  conversion  of  ethylene  and  the 
yield  of  ethylene  oxide  at  the  outlet  of  the  reactor,  re- 
spectively. 

The  optimization  problem  is  to  maximize  the  total  amount 
of  yield  of  ethylene  oxide  over  a  fixed  total  reaction 
time,  tj,  by  choosing  the  optimal  temperature-time  policy. 
We  assume  that  the  temperature  dependence  of  the  rate  con- 
stants and  adsorption  equilirbium  constant  are  all  of  the 
Arrhennius  form: 

K±  =  Ai  exp(-Ei/RgT)    i  =  1,  2,  q,  D     (5-6) 

Since  K  is  a  strictly  monotonic  function  of  temperature, 

we  can  use  KD  as  the  control  variable  to  replace  T.   The 

constants  K, ,  K„  and  K   can  be  expressed  in  terms  of  K^: 
J-    t.  q  D 
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K.  =  a.K 


where 


J  D 


a.  =  A./V 


J" 


1,2, q 


(5-7) 


and 


Pj  =  Ej/ED 


Thus,  our  optimization  problem  is 


r  t. 


J  =  max  • 
KD(t) 


Ydt  \   =   max 
KD(t) 


K. 


lJ0 


Kl  +  K2 


Xdtf    (5-8) 


subject  to  two  local  restrictions,  Eqs.  (5-3)  and  (5-4). 
This  is  a  fundamental  problem  in  calculus  of  variations. 
The  derivation  of  the  necessary  conditions  for  J*  to  be 
maximum  is  presented  in  the  next  section. 
§5-2  Necessary  Condition 

Equation  (5-3)  indicates  that 


where 


KD  =  KD(h,  h) 


h  =  dt 


(5-9) 


Since  1^,  K2  and  K  are  all  functions  of  K  ,  Eqs.  (5-4) 
and  (5-5)  give 


Y  =  Y(KD,  h)  =  Y(h,  h) 


(5-10) 


Thus,  the  optimization  problem  of  Eq.  (5-8)  is  an  elemen- 
tary problem  in  the  calculus  of  variations,  and  the  necessary 
condition  for  the  optimization  is  an  Euler-Lagrange   equa- 
tion (Beveridge  and  Schechter,  1970): 


where 


Y  -  h  Y; 


8Y 


(5-11) 


3h 


h,t 


and  C  is  an  integration  constant. 
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In  order  to  use  Eq.  (5-11),  hY^  has  to  be  obtained 
first.   Differentiating  Eq.  (5-7)  w.r.t.  h  yields 


3K, 


3h 


rdK.  > 

l 

f9V 

*>    3h' 

= 

L    KD    J 

r    i  | 

k    h      J 

i  =  1,  2,  q    (5-12) 


Multiplying  h  gives 


3K. 

l 

3n 


rp.K. 


K, 


_h 
m 


=  p.K. 


(5-13) 


Differentiating  Eq.  (5-4)  w.r.t.  h,  multiplying  n  and 
using  Eq.  (5-13)  yields 

PqV  +  V*  Xh)  +  T^X  ^  Xh}  =  P3T(K1  +  K2)h    (5"14) 
where     p3  =  (p1K1  +  p2K2)/(K1  +  K2 ) 

Rearranging  Eq.  (5-14)  and  using  Eq.  (5-4)  for  t (K   +  K  )h 
we  have 


.        (P3~  P  )K   X  -  p   An(l-X) 

h  JO  =  — = 3 — 9 2 

h 

q    1-X 


(5-15) 


Differentiating  Eq.  (5-5)  w.r.t.  h,  multiplying  h  and 
using  Eq.  (5-13)  we  obtain 
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K  K         .  K 

h  Yh =  ?!   k— nqx  +  K7inqh  xh  -  p3  k— tk7x 

(5-16) 
Substituting  Eq.  (5-16)  into  (5-11)  with  the  aid  of 
Eq.  (5-15)  gives 


K 
K~ 


K 

P Tiq   K   731   {(V  P3,KqX  +  P3ta<l-«  +  V  1^J 

q    1-X 

(1-Pl  +  p_)xl  =  C    (5-17) 


Equation  (5-17)  is  the  necessary  condition  for  J  to  be 
maximum.   The  integration  constant  C  can  be  determined  by 
a  boundary  condition: 

h  =  hf   at   t  =  tf  (5-18) 

if  h^  is  specified,  or  by  a  natural  boundary  condition 
(Denn,  1969) : 

Y^  =  0     at     t  =  tf  (5-19) 

if  hf  is  unspecified. 

The  necessary  condition  can  be  simplified  if  the  acti- 
vation energies  for  both  reactions  are  equal  (E.  =  E~) , 
and  the  adsorption  equilibrium  constant  K   is  temperature 

q 

independent  (E  =0).   In  this  case,  Eq.  (5-17)  gives 


Al 

J-p..A  X  +  p,£n(l-X)  +  (A   +  t^;)x}-  =  C 


(A1  +  A2)(Aq  +  TTx>  l   ^  q    ^ 

(5-20) 
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where  A^   is  the  preexponential  factor.   Since  A  ,  A  ,  Aq, 
P1  and  C  are  all  constants,  Eq.  (5-20)  indicates  that 

X  =  constant  (5-21) 

In  general  cases,  however,  the  conversion  is  not  constant. 
The  optimal  control  policies  have  to  be  obtained  by  numeri- 
cally solving  the  necessary  condition,  Eq.  (5-17)  ,  along 
with  the  local  constraints,  Eqs.  (5-3)  and  (5-4). 
§5-3  Optimal  Continuous  Control  Policies 

The  predictor-corrector  method  (Hornbeck,  1975)  initiated 
by  fourth  order  Runge  Kutta  method  was  used  to  solve  the 
ordinary  differential  Eq.  (5-3)  while  Newton-Raphson  itera- 
tion (Carnahan  et  al.,  1969)  was  used  for  the  nonlinear 
Eqs.  (5-4)  and  (5-17).   The  kinetic  data  obtained  in  Section 
§3-2-3  and  the  same  parameters  used  in  the  experiment  were 
used  for  obtaining  the  optimal  control  policies.   Table  5-1 
lists  these  parameters  along  with  the  deactivation  kinetics 
used.   The  boundary  condition  of  Eq.  (5-18)  was  used  to 
obtain  the  integration  constant  in  Eq.  (5-17).   This  split 
boundary  condition  problem  was  solved  by  "shooting  method" 
(Hornbeck,  1975).   The  results  are  shown  in  Figs.  5-1  and 
5-2. 

It  can  be  seen  from  Fig.  5-1  that  the  optimal  control 
policy  is  one  of  increasing  conversion  and  yield.   However, 
with  different  deactivation  kinetics,  the  optimal  policy 
becomes  one  of  decreasing  conversion  and  yield,  as  shown 
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Table  5-1   Parameters  Used  for  the  Optimal  Continuous  Con- 
trol Policies 

Kl  =  klkq  =  7'028  x  10~2  exp(-5451/T) (mol/s • g-cat - atm) 
K2  =  k2kq  =  5-900  x  10-1  exp(-6632/T) (mol/s -g-cat-atm) 
Kq  =  PAokq  =  1'165  X  10_1  exp(1647/T)  (-) 


x   ~   p~  =  1.134  x  10        ( g-cat • atm -s /mo 1) 
o 

hf  =  0.4   (-) 
tf  =  1000  (hr) 

KD  =  4.0  x  102  exp(-6000/T)  (1/hr)  in  Fig.  5-1 

=  4.0  x  103  exp(-7000/T)  (1/hr)  in  Fig.  5-2 

El 
?!=£-=  0.9085  in  Fig.  5-1 

=  0.7787  in  Fig.  5-2 

E2 
P2  =  e~  =  1-105  in  Fig.  5-1 

=  0.9474  in  Fig.  5-2 

m  =  l.o 


El/E2  =  °-822 
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in  Fig.  5-2.   This  decrease  in  conversion  and  yield  is 
due  to  the  large  activation  energy  for  deactivation  which 
limits  the  increase  in  temperature. 

The  valuesof  the  performance  index  for  constant-conver- 
sion and  constant-yield  policies  have  been  calculated  for 
comparison.   The  results  indicate  that  the  control  policies 
shown  in  Figs.  5-1  and  5-2  are  indeed  the  optimal  policies. 
However,  the  percent  difference  between  the  optimal  policy 
and  the  constant-conversion  policy  is  less  than  0.1%  in 
terms  of  the  performance  index.   As  mentioned  earlier,  the 
constant-conversion  policy  is  the  optimal  policy  when  E  =  E„ 
and  E  =  0.   For  the  model  reaction  studied,  E1    is  close  to 
E2  (E1/E2  =  0.822)  and  E   is  small  (compared  with  E   and 
E„) .   This  explains  why  the  constant-conversion  policy  is 
almost  as  good  as  the  optimal  policies  shown  in  Figs.  5-1 
and  5-2. 
§5-4   Concluding  Remarks 

We  have  derived  the  necessary  condition  for  the  optimal 
continuous  control  for  a  parallel  reaction  network  with 
Langmiur-Hinshelwood  rate  expressions  by  direct  substitution. 
The  optimal  control  policies  have  been  obtained  by  direct  nu- 
merical integration  of  Eq.  (5-3)  with  the  aid  of  Eqs.  (5-4) 
and  (5-17) .   Although  iterations  are  required  for  solving 
Eqs.  (5-4)  and  (5-17),  it  takes  no  more  than  two  iterations 
for  each  step  in  the  t-direction  because  X  and  K   are  smooth 
functions  of  t  and  the  values  in  the  previous  step  provide 
very  good  initial  guesses  for  the  current  step.   In  fact,  it 
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took  less  than  3  seconds  CPU  time  of  IBM  360/370  computer 
to  generate  all  the  data  shown  in  Fig.  5-1,  including  5 
iterations  for  shooting  the  boundary  condition,  h  =  h,  at 
t  =  tf .   This  efficiency  can  never  be  achieved  if  the 
Pontryagin  maximum  principle  (Pontryagin  et  al.,  1962) 
and  the  Hamiltonian  are  used  (Ogunye  and  Ray,  1971) . 


CHAPTER  6 
OPTIMAL  CONTROL  AND  DESIGN 
We  address  in  this  chapter  the  problem  of  how  a  fixed- 
bed  with  catalyst  deactivation  should  be  designed  to 
obtain  the  best  possible  performance  when  it  is  controlled 
optimally.   In  particular,  we  ask  what  the  size  of  the 
reactor  should  be  that  maximizes  a  certain  performance  index 
when  the  temperature  is  also  manipulated  to  maximize  the  per- 
formance index.   This  is  a  design  practice  in  which  process 
control  is  taken  into  consideration.   In  order  to  clearly  de- 
monstrate the  interrelationship  between  design  and  control 
and  the  advantage  that  can  be  gained  by  utilizing  this 
interrelationship,  we  assume  for  the  time  being  that  the 
detailed  knowledge  of  deactivation  is  available. 
§6-1   An  Optimization  Problem  of  Piecewise  Control  and  Design 

Consider  an  optimization  problem  of  a  single  irre- 
versible reaction  taking  place  in  an  adiabatic  reactor 
given  by 

*  ,  rt 


J   =   max       J  = 


max 


T'  Tin(t)         T'  Tin(t)   Ttf 


f   Xdt 


(6-1) 

subject  to  the  system  equations  of  Eqs.  (2-1)  and  (2-2) 

with  (mCp)c  =  0  for  the  adiabatic  reactor  under  consideration 
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and   subject  to     T^n  <.    (Tj_n)max  with  a  fixed  t_ .   The 
problem  is  to  maximize  the  performance  index  J,  which  is 
the  conversion  (X)  integrated  over  the  catalyst  life  tf 
per  unit  reactor  volume   and  unit  on-stream  time,  by  proper- 
ly choosing  x  and  Tin(t).   This  two-parameter  optimization 
problem  can  be  solved  one  at  a  time  by  first  searching 
for  the  maximum  for  a  given  x  and  then  finding  the  value  of 
x  that  yields  the  maximum.   For  a  fixed  x,  it  is  a  Bolza's 
problem  (Bliss,  1961)  and  the  extremum  condition  obtained 
by  Chou  et  al.  (1967)  and  Ogunye  and  Ray  (1971)  is  to  main- 
tain the  conversion  at  a  constant  level  until  the  maximum 
allowed  temperature  is  reached  and  then  to  remain  at  the 
maximum  if  only  one  rate  constant  is  involved.   For  the  de- 
activation problem  under  consideration  the  policy  is  not 
directly  applicable  since  the  rate  constant  cannot  be  se- 
parated out.   However,  the  constant  conversion  policy  is 
still  a  suboptimum  which  is  often  used  in  practical  opera- 
tion of  chemical  plants.   Thus,  we  apply  the  constant 
conversion  policy  for  the  problem  considered.   Since 
the  feedback  control  policy  obtained  earlier  is  for  main- 
taining a  constant  conversion  level,  the  problem  is  simply 
that  of  finding  an  initial  value  of  T.   (or  the  desired 
conversion  level)  if  we  use  the  control  policy,  for  T 


in 


can  be  set  at  (Tin)max  once  (Tin)max  is  reached 


75 

The  optimization  problem  was  solved  for  the  reaction 

system  in  Tables  3-1  and  3-2  using  the  control  algorithms 

of  Eqs.  (2-16)  through  (2-18)  with  the  constraints  of  tf  = 

10   sec.  and  (T.  )     =  720  K.   The  typical  behavior  of 
in  max 

T.  (t)  and  the  outlet  conversion  is  shown  in  Fig.  6-1  for 

the  bandwidth  AX  of  5%  and  the  initial  inlet  temperature 

of  691  K  (or  conversion  of  87.2%).   In  accordance  with 

the  constraint  on  t_  and  the  extremum  condition,  T.   stays 

f  in    ^ 

at  the  maximum  allowed  temperature  once  it  reaches  the  maxi- 
mum and  the  conversion  decreases  with  time,  eventually 
reaching  a  conversion  of  54.5%  at  t  =  tf,  at  which  time 
the  catalyst  is  regenerated.   As  indicated  earlier,  the 
optimization  problem  for  a  chosen  x  is  that  of  finding  an 
initial  conversion  (or  initial  inlet  temperature)  that  yields 
a  maximum.   The  solutions  obtained  for  various  t  are  shown 
in  Fig.  6-2  for  a  bandwidth  of  0.05.   A  few  observations 
can  be  made  from  the  figure.   The  results  show  that  for  a 
given  x,  the  extremum  condition  does  yield  a  maximum. 
Furthermore,  the  value  of  the  performance  index  is  much 
more  sensitive  to  a  choice  of  the  initial  inlet  temperature 
(or  the  initial  conversion)  at  relatively  lower  values  of 
the  inlet  temperature  than  at  higher  values.   The  locus 
of  the  maximum  values  of  J  for  various  x  is  given  in 
Fig.  6-2  as  the  dotted  line.   As  apparent  from  the  figure, 
there  exists  a  maximum  at  around  x  =  22  sec.   The  perfor- 
mance index  normalized  with  respect  to  the  index  corresponding 
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Figure  6-2.   Performance  Index  as  a  Function  of 
Reactor  Size  and  Initial  Inlet 
Temperature.   The  dotted  line  shows 
the  locus  of  J*. 
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to   t   =   22  sec,  when  the  reactor  is  operated  optimally, 
is  shown  in  Fig.  6-3  as  a  function  of  the  reactor  size  t. 
It  is  seen  that  the  performance  index  [Eq.  (6-1)]  doubles 
(t  of  8  sec.  vs.  t  of  22  sec.)  when  the  reactor  is  sized 
optimally.   This  result  confirms  the  intuitive  reasoning 
that  the  best  possible  reactor  performance  obtainable 
by  process  control  is  inherently  dictated  by  process  design, 
It  also  confirms  the  fact  that  a  significant  improvement 
can  be  made  by  combining  process  design  with  process  con- 
trol. 

§6-2   An  Optimization  Problem  of  Continuous  Control  and 
Design  " 

In  this  section,  we  give  an  example  of  design  practice 
in  which  the  optimal  continuous  control  is  taken  into 
consideration.   Here  we  consider  an  isothermal  fixed-bed 
where  a  first  order  reaction  subject  to  first  order  deacti- 
vation takes  place.   The  system  equations  are  the  same  as 
those  in  Table  4-1,  but  the  temperature  is  controlled  con- 
tinuously.  The  performance  index  is  given  by 


j*  _  max         max 


-  *f 


T,T(t)        T,T(t)    Tt4 


Xdt  -  y 
0 


(6-2) 

where  y  is  the  cost  incurred  due  to  the  reactor  shut-down 
for  regenerating  or  replacing  the  spent  catalyst.   The 
constant  conversion  policy  which  is  indeed  the  optimal 
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control  policy  for  this  reaction  system  (Crowe,  1970) 

was  applied  to  solve  the  optimization  problem  of  Eq.  (6-2) 

with  the  following  rate  constants  and  restrictions 

kD  =  1.0  x  104  exp(-8000/T)      (1/hr) 

kx  =  1.0  x  106  exp(-10000/T)     (1/S) 

473  K  <      T  <_      673  K 

tf  =  2000  (hr) 

y  =  220  (hr) 

The  solutions  obtained  are  plotted  in  Fig.  6-4.   The 

behavior  of  the  performance  index  is  similar  to  that 

shown  in  Fig.  6-2,  i.e.  for  a  given  reactor  size  (x)  ,  J 

depends  on  the  initial  temperature  chosen  (T  )  and  does 

o 

yield  a  maximum.   The  locus  of  the  maximum  value  of  J  for 
various  of  x  is  shown  in  Fig.  6-4  as  the  dotted  line. 
Again,  we  can  see  from  the  figure  there  exists  an  optimal 
reactor  size,  i.e.  x  =  13  sec.   Further  observations  of  the 
figure  indicate  that  a  reactor  of  x  =  15  second  gives 
better  performance  than  a  reactor  of  x  =  10  second  does  if 
both  are  controlled  optimally   (T   =  595.9  K  for  x  =  15 
sec.  and  Tq  =  605.4Kfor  x  =  10  sec).   However,  the  latter 
performs  better  than  the  former  does  if  same  initial  tem- 
peratures are  chosen,  say,  620  K.   This  indicates  that 
both  the  deactivation  of  catalyst  and  the  corresponding 
process  control  should  be  considered  during  the  design  phase 
to  choose  an  optimal  reactor  size. 
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Figure  6-4. 


Performance  Index  "as  a  Function  of 
Reactor  Size  and  Initial  Temperature* 
The  dotted  line  shows  the  locus  of  J 
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The normalized  optimal  performance  index  is  plotted  in 
Fig.  6-5  as  a  function  of  reactor  size.   It  can  be  seen 
from  the  figure  that  the  performance  index  can  be  improved 
significantly  if  the  optimal  size  is  chosen. 
§6-3   Concluding  Remarks 

We  have  so  far  assumed  that  detailed  knowledge  of  ca- 
talyst deactivation  is  available.   The  same  optimal  con- 
trol and  design  can  still  be  carried  out  if  some  knowledge 
of  deactivation  is  available,  however  uncertain  it  may  be, 
since  what  is  sought  in  such  a  case  would  be  approximate 
values  of  the  optimal  size  and  inlet  temperature  and  some 
directions  as  to  what  to  choose.   In  the  example  considered 
in  Section  §6-1,  for  example,  one  would  choose  a  relative- 
ly high  initial  conversion  (or  a  high  inlet  temperature) 
rather  than  a  low  initial  conversion  in  view  of  the  sen- 
sitivity discussed  earlier. 
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CHAPTER  7 
CONCLUSIONS 
A  measure  of  the  extent  of  deactivation  can  be  deter- 
mined from  the  on-line  measurements  of  inlet  temperature 
and  outlet  concentration  along  with  the  intrinsic  rate  of 
a  single  reaction  for  an  isothermal  or  an  adiabatic  fixed- 
bed  reactor.   Additional  measurements  of  coolant  tempera- 
ture are  required  for  a  nonadiabatic  fixed-bed  reactor. 
These  results  can  be  extended  to  multiple  reactions  in  a 
rather  straightforward  manner.   The  measure  of  deactivation 
is  determined  in  such  a  way  that  a  piecewise  feedback  con- 
trol scheme  follows  directly  from  it.   This  feedback  con- 
trol maintains  the  conversion  within  a  band  about  a  desired 
level  of  conversion.   These  results  enable  one  to  carry  out 
feedback  control  of  the  reactor  in  any  desired  manner  with- 
out any  knowledge  of  catalyst  deactivation.   Further,  the 
results  are  applicable  to  any  type  of  deactivation;  whether 
the  deactivation  is  dependent  or  independent  of  reactant 
concentration  and  whether  the  deactivation  kinetics  are 
separable  or  nonseparable  from  main  reaction  kinetics. 

The  feedback  control  policy  considered  is  a  piecewise 
rather  than  continuous  control.  The  piecewise  control  ne- 
cessitated by  the  pseudo  steady-state  nature  of  deactivation 
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with  respect  to  the  reactor  response  time  allows  us  to  avoid 
one  problem  in  any  feedback  control  involving  dead  time  due 
to  sampling  and  another  one  associated  with  reactor  tran- 
sients.  Inherent  in  the  pseudo  steady-state  assumption  is 
the  fact  that  the  rate  of  deactivation  is  much  slower  than 
the  rate  of  reaction.   In  terms  of  time,  the  deactivation  is 
of  the  order  of  months  for  the  majority  of  reactions  carried 
out  in  a  fixed  bed  whereas  the  transients  reach  the  steady 
state  in  a  few  residence-times  which  is  of  the  order  of 
seconds.   Therefore,  the  problem  associated  with  "wrong  way" 
transients  [e.g.,  Crider  and  Foss  (1966)]  can  be  avoided 
since  by  the  time  a  control  action  is  taken,  the  reactor 
reaches  a  new,  stable  steady  state.   Further,  the  effect  of 
dead  time,  which  is  of  the  order  of  minutes,  on  deactivation 
can  be  neglected  since  the  extent  of  deactivation  changes 
little  in  minutes. 

The  feedback  control  policy  has  been  applied  to  a 
simulated  adiabatic  reactor  where  a  single  reaction  affected 
by  diffusion  and  poisoning  take  place,  and  to  a  laboratory 
isothermal  reactor  where  ethylene  oxidation  reactions  affect- 
ed by  sintering  take  place.   The  results  show  that  one  can 
indeed  control  the  conversion  or  yield  in  any  desired  manner 
without  any  a  priori  knowledge  of  the  deactivation  kinetics, 
and  that  the  aforementioned  problems  of  wrong  way  transients 
and  dead  time  do  not  occur  during  normal  operation. 
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The  problem  of  finding  the  optimal  piecewise  con- 
trol policy  can  be  formulated  as  a  multivariable  optimi- 
zation problem.   For  a  single  irreversible  reaction  affected 
by  independent  deactivation,  the  optimal  piecewise  control 
policy  is  one  of  increasing  conversion.   However,  the  well- 
known  constant  conversion  policy  used  in  the  continuous  con- 
trol system  can  be  applied  to  the  piecewise  control  sys- 
tem with  little  loss  in  reactor  performance  if  the  number  of 
piecewise  control  actions  is  large.   The  solution  method 
developed  is  applicable  to  any  reaction  taking  place  in 
any  reactor. 

The  necessary  condition  for  optimal  control  in  complex 
reaction  systems  (such  as  ethylene  oxidation  reactions) 
can  be  derived  by  direct  substitution.   The  optimal  control 
policy  is  either  an  increasing  conversion  policy  or  a  de- 
creasing conversion  policy,  depending  upon  the  ratio  of 
activation  energy  of  the  main  reaction  to  that  of  the 
deactivation  reaction.   The  direct  substitution  method  is 
shown  to  be  more  efficient  than  the  commonly  used  method 
that  relies  on  the  Pontryagin's  minimum  principle  through 
a  Hamiltonian  formulation. 

The  problem  of  combining  process  control  with  process 
design,  with  due  consideration  of  catalyst  deactivation 
for  both,  can  be  solved  in  a  straightforward  manner  using 
the  feedback  control  scheme  and  the  approach  of  reactor  point 
effectiveness . 
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However,  detailed  knowledge  of  deactivation  has  to  be  known 
to  solve  this  problem.   It  is  shown  that  the  choice  of  re- 
actor size  can  substantially  affect  the  reactor  performance. 
The  solution  yields  the  optimal  reactor  size  and  an  optimal 
way  of  manipulating  the  inlet  temperature  based  on  on-line 
measurements  for  those  reactors  for  which  detailed  know- 
ledge of  deactivation  is  available.   Even  when  there  is 
some  uncertainty  about  the  knowledge,  the  procedure  yields 
approximate  ranges  of  the  optimal  values  and  also  some  di- 
rections as  to  what  is  to  be  chosen.   These  results  enable 
one  to  fully  account  for  catalyst  deactivation  for  a  substan- 
tial improvement  in  the  reactor  performance. 
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